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ABSTRACT 

With  long  pulse-duration  (4  nsec),  moderate  intensity  (=1013  W/cm2)  Nd-laser  il¬ 
lumination  of  planar  targets,  evidence  is  found  in  images  of  x-ray,  3/2  w0  and  2&>0  emis¬ 
sion  for  filamentation  of  the  incident  beam.  Calculations  indicate  that  this  may  be  due  to 
thermal  self-focusing  of  hot  spots  in  the  laser.  X-ray  images  show  spots  and  spicules,  jets 
of  material  flow  which  appear  to  follow  the  fluid  flow  from  the  target.  These  may  be  at¬ 
tributed  to  the  presence  of  microscopic  impurities  on  our  laser-irradiated  foil  targets. 

NOTE:  This  manuscript  has  been  revised  in  proof  to  include  information,  as  indicat¬ 
ed  in  Sec.  IV,  that  was  not  available  at  the  time  of  the  oral  conference  presentation.  This 
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information  is  also  contained  in  a  published  erratum. 
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I.  INTRODUCTION 

High  compression  of  laser-fusion  pellets 
requires  application  of  a  uniform  pressure  at  the 
pellet  surface.1  While  driving-pressure  symmetry 
can  be  enhanced  due  to  lateral  energy  flow  during 
transport  between  the  absorption  region  and  the 
pellet  surface,1  2  it  is  certainly  desirable  for  the 
laser  illumination  in  the  absorption  region  to  be  as 
uniform  as  possible.  As  presently  envisioned, 
multimillimeter-diameter  pellets  will  be  irradiated 
by  multinanosecond  pulses  in  laser  fusion  reactors. 
In  the  large  plasma  atmospheres  that  will  result, 
thermal3-6  or  ponderomotive3  7-10  effects  may  lead 
to  self-focusing  or  filamentation  of  the  incident 
laser  beam,  thereby  seriously  degrading  illumina¬ 
tion  uniformity.  In  addition,  other  effects  that  are 
deleterious  to  laser  fusion  (e.g.,  fast  electron  gen¬ 
eration)  may  be  driven  locally  by  the  higher  laser 
intensity  present  in  the  filaments. 

To  simulate  the  large  plasma  atmosphere  of  a 
reactor  pellet,  50-300  Joule  pulses  from  the  Pharos 
II  Nd-laser11  (A  =  1.054  gm,  tl  =  3-4  nsec)  are 
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focused  to  large  spot  diameter  (0.6-1. 0  mm)  on 
planar  targets  of  polystyrene  (CH)  or  aluminum. 
Large  underdense  plasmas  (lengths  near  1  mm, 
with  density  gradient  scalelengths  near  ne  =  0.1  nc  in 
excess  of  100  /xm12)  are  thereby  created.  Images  of 
second  harmonic  (2tu0)  and  three-halves  subhar¬ 
monic  (3/2  ojq)  emissions  provide  evidence  for 
filamentation,  these  are  discussed  in  Sec.  III.  X- 
ray  images  show  spots13  and  spicules,14  jets  of 
material  flowing  from  the  target  surface.  These 
may  be  attributed  to  the  presence  of  microscopic, 
higher-Z  impurities  on  the  low-Z  foil  targets. 


II.  EXPERIMENTAL  DESCRIPTION 

As  stated  in  the  introduction,  planar  targets  of 
CH  or  At  are  illuminated  with  50-300  J  of  Nd-laser 
(A  =  1.054  /xm,  tl  =  3-4  nsec)  radiation.  In  the 
quasi-near-field  of  an  f/7  input  lens,  spot  diameters 
(50%  energy  content)  between  0.6  and  1.0  mm  are 
irradiated  with  incident  intensities  /0  =  0.1-1  x 
1013  W/cm2.  Targets  are  tilted  by  6°  relative  to  the 
incident  beam  direction  so  as  to  deflect  specular 
reflections  from  the  incident  lens  cone.  Incident 
beam  polarizations  are  either  circularly  or  p- 
poiarized;  since  no  major  differences  are  found  in 
the  results  for  these  two  polarizations,  no  further 
mention  of  polarization  state  will  be  made. 
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The  experimental  arrangement  is  schemati¬ 
cally  shown  in  Fig.  1.  Visible  emission  is  collected 
parallel  to  the  target  surface  by  an  f/2  lens  and 
imaged  through  narrow  band  filters  onto  Polaroid 
positive-negative  film  (for  2to0)  or  Kodak  type  1-N 
spectroscopic  plates  (for  3/2  to0).  These  images, 
obtained  with  the  "2-D  camera"  of  Fig.  1,  yield 
time-integrated  information  with  two  degrees  of 
spatial  resolution.  A  Wollaston  prism  in  the  imag¬ 
ing  path  causes  slight  degradation  in  focal  sharp¬ 
ness,  but  allows  simultaneous,  independent  imaging 
in  two  polarizations  of  emitted  light.  This  enhances 
contrast  of  polarized  harmonic  emission  over  unpo¬ 
larized  continuum  emission.  Time-resolved  infor¬ 
mation  can  also  be  obtained  by  imaging  a  fraction 
of  the  emission  onto  the  slit  of  a  streak  camera, 
though  only  one  spatial  dimension  is  resolved  in 
these  streak  images. 


filters,  so  that  images  are  of  photons  with  energy 
hi’  >  1  keV. 


III.  EVIDENCE  FOR  FILAMENTATION 

In  the  present  work,  experimental  evidence  is 
found  for  four  signatures  of  filamentation:  (1) 
localization  on  a  finer  spatial  scale  than  incident 
beam  nonuniformities,  (2)  locally  lower  density  in 
these  regions,  (3)  locally  higher  laser  intensity  in 
these  regions,  with  an  increase  in  intensity  modula¬ 
tion  over  that  present  in  the  incident  beam,  and  (4) 
the  appearance  of  these  effects  only  at  higher 
incident  intensity.  Previous15  16  and  concurrent17  18 
investigations  provide  evidence  for  some  of  these 
properties;  the  simultaneous  inference  of  all  four  in 
this  work  allows  a  stronger  case  to  be  made  for  the 
presence  of  filamentation. 


EXPERIMENTAL  APPARATUS 


A.  Uniformity  of  the  Incident  Beam 


2  D  CAMERA 


Fig.  I  —  Experimental  schematic,  showing  imaging  systems 
for  visible  and  x-ray  emission.  All  visible  imaging  was 
done  from  the  side. 


Time-integrated  images  of  x-ray  emission 
from  these  targets  are  obtained  from  a  range  of 
angles,  as  shown  in  Fig.  1.  Both  frontal  and  rear 
pinhole  cameras  detect  emission  from  the  laser- 
irradiated  front  side  of  the  target,  though  the  rear 
camera  looks  only  at  emission  which  is  transmilt  j 
through  the  target  material.  All  cameras  are 
equipped  with  nominally  18  ftm  thick  beryllium 


We  must  characterize  the  nonuniformity  of 
the  focused  incident  beam  in  order  to  recognize 
changes  in  the  intensity  distribution  which  are 
caused  by  interaction  with  the  plasma.  In  Fig.  2B  is 
shown  an  example  of  one  isointensity  contour  from 
a  series  of  contours  which  is  obtained  on  bismuth- 
coated  glass  slides  in  an  equivalent  focal  plane  on 
every  target  shot.  Large  scale  structure  in  the  focal 
distribution  is  seen  in  the  form  of  linear  fringes 
with  100-200  jxm  spacing;  inspection  of  the  entire 
series  of  isointensity  images  suggests  a  peak-to- 
trough  intensity  modulation  level  near  2:1  for  this 
structure.  Interestingly,  this  pattern  is  also 
observed  in  x-ray  images  as  obtained  with  the  rear 
pinhole  camera,  as  shown  in  Fig.  2A. 

The  finest  structure  observed  in  laser  focal 
images  has  50-80  spatial  scale  and  peak  intensi¬ 
ties  in  these  hot  spots  may  be  up  to  twice  the  aver¬ 
age  intensity;  while  even  finer  structure  may  exist, 
it  must  have  significantly  lower  modulation  level. 
Finer  scale  structure  observed  in  x-ray  images,  as 
best  illustrated  by  Fig.  2C,  was  previously 
believed13  to  provide  evidence  for  more  localized 
regions  of  higher  laser  intensity;  however,  we  now 
have  evidence,  as  discussed  in  Sec.  IV,  for  the 
presence  of  microscopic,  high-Z  impurities  in  our 
low-Z  target  foils.  While  we  cannot  state  with  cer¬ 
tainty  that  all  of  the  fine  scale  structure  observed  in 
x-ray  images  such  as  Fig.  2C  is  due  to  these  impuri¬ 
ties,  their  presence  precludes  our  use  of  such  x-ray 
images  as  evidence  for  laser  intensity  modulation. 
The  most  direct  evidence  for  filamentation  is  based 
upon  the  visible  observations  discussed  in  Sec.  Ill 
B. 
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Fig.  2  —  (A)  Rear-side  x-ray  image  for  12.5  (im  Al  foil  target  at  /0  -  2.9  x  lo12  W/cm2;  (B) 
laser  focal  isointensity  contour  from  same  shot  as  (A);  (C)  Rear  side  x-ray  image  for  7  /am  Al 
foil  target  at  /0  =  3.4  x  1012  W/cm2. 


B.  Visible  Observations 

As  the  incident  intensity  is  increased,  time- 
integrated  images  of  emitted  light  at  the  second 
harmonic  (2<«0)  and  three-halves  subharmonic  (3/2 
a/0)  reveal  that  emission  first  occurs  from 
filamentary-appearing  regions,  with  complete  sur¬ 
faces  of  emission  (presumably  the  ncl 4  surface  for 
3/2  to0  and  the  nc  surface  for  2oj0)  appearing  only 
at  higher  incident  intensity.  Near  threshold  for 
detection  of  3/2  <o0  radiation  (average  intensity 
/0  =  5.6  x  !012  W/cm2),  filamentary  structures  as 
in  Fig.  3A  are  seen  in  both  emitted  light  polariza¬ 
tions.  These  filaments  have  transverse  dimensions 
of  5-15  gm  and  lengths  greater  than  200  gm.  As 
/0  is  increased  by  a  factor  of  2,  3/2  w0  emission 
from  the  entire  ncl 4  surface  is  detected;  this  emis¬ 
sion  appears  first  on  the  downstream  (laser)  side  of 
the  filaments,  but  is  quite  broad  in  extent  along  the 
laser  direction.  Two  lobes  of  emission  are  observed 
in  images  of  3/2  «o0  with  polarization  corresponding 
to  the  laser  (Fig.  3B);  curiously,  a  different  image 
is  obtained  in  the  orthogonal  polarization  (Fig.  30. 

With  /0  >  2.8  x  1012  W/cm2,  filamentary 
structures  are  also  found  in  time-integrated  images 
of  2  (u0  emission,  as  seen  in  Fig.  4.  The  large 
region  of  uniform  exposure  is  visible  continuum,  as 
verified  by  substitution  of  filters  which  exclude 
2  ai0;  only  the  filaments  along  the  laser  axis  are 
unique  to  2  w0.  The  filaments  have  lengths  along 
the  laser  axis  of  up  to  125  gm,  transverse  dimen¬ 
sions  of  5-15  gim,  and  are  observed  only  in  images 
with  polarization  corresponding  to  the  incident 
laser.  With  a  factor  of  4  increase  in  /0,  parts  of  a 
well-defined  surface  of  2  <a o  emission  appear  at  the 
downstream  end  of  the  filaments. 

Inference  of  Local  Intensity  Modulation.  From 
the  transverse  modulation  of  the  3/2  o>0  emission, 


we  infer  the  presence  of  localized  regions  of  higher 
laser  intensity,  with  an  intensity  modulation  level  at 
least  as  large  as  that  present  in  the  incident  h-'am. 
As  a  complete  surface  of  3/2  &>0  emission  is 
observed  at  twice  the  intensity  at  which  the 
filaments  are  first  seen,  one  might  infer  that  the 
intensity  in  the  filaments  is  twice  that  in  the  sur¬ 
rounding  beam.  However,  the  conversion  efficiency 
of  light  from  w0  to  harmonics  is  a  function  of  den¬ 
sity  gradient  scalelength  L„  =  (nf'  dnj dx)~ ',  as 
well  as  of  laser  intensity.  The  efficiency  of  conver¬ 
sion  from  o»0  to  3/2  o>0  increases  monotonically  with 
density  gradient  scalelength19-21  (at  least  until 
scalelengths  become  so  long  that  inverse 
bremsstrahlung  reduces  the  amount  of  light  reach¬ 
ing  the  nc/4  region21).  As  one  would  expect  local¬ 
ized  structures  to  have  shorter  scalelengths  than  the 
rest  of  the  blowoff,  the  3/2  w0  production  would  be 
reduced,  if  anything,  for  a  given  incident  intensity 
in  the  observed  channels.  We  may  then  estimate  a 
lower  bound  for  the  intensity  modulation  at  nJ4  by 
assuming  that  Ln  in  the  filaments  and  in  the 
remainder  of  the  blowoff  are  equal;  this  lower 
bound  is  the  2:1  intensity  modulation  level  men¬ 
tioned  above.  The  lower  bound  is  the  same  as  the 
intensity  modulation  level  in  the  incident  beam,  but 
there  is  now  a  finer  spatial  scale. 

Similar  analysis  of  the  2  w0  images  yields  an 
intensity  modulation  at  «,  of  4:1,  as  parts  of  a  sur¬ 
face  are  seen  at  4  times  the  incident  intensity  at 
which  the  filaments  are  seen.  Theoretically,  the 
dependence  of  2  to0  generation  on  L„  is  nonmono¬ 
tonic,  however,  maximizing  at  density  gradient 
scalelengths  of  several  wavelengths.22  24  Since  it  is 
possible,  therefore,  to  have  enhanced  2  w0  genera¬ 
tion  due  to  shorter  L„,  it  is  not  clear  whether  the 
inferred  intensity  modulation  is  an  overestimate  or 
underestimate.  A  4:1  modulation  would  certainly 
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Fig  3  —  3/2  o)0  images  for  CH  targets:  (A)  central  region  of  I  mm  laser  spot  on  40  ^m  thick 
foil  at  /0  =  6.1  x  1012  W/cm2;  (B)  42.7  fim  thick,  500  /am  diameter  disk  at  /0  —  l  lx  1013 
W/cm2,  with  emitted  light  polarization  corresponding  to  incident  laser;  and  (C)  same  shot  as 
(B),  but  with  orthogonal  polarization.  Laser  is  incident  from  left  in  all  images. 


Fig.  4  —  2o)0  image  for  7  /im  CH  target  at  l0  =  6.1  x  1012 
W/cm2.  Laser  is  incident  from  left,  and  dark  vertical  edge 
at  right  is  target  surface.  General  exposure  due  to  emitted 
continuum  is  seen  in  addition  to  filaments  of  emitted  2<u0. 


represent  an  increase  over  the  modulation  present 
on  the  incident  beam. 

Mechanisms  That  Can  Spuriously  Cause  Fila¬ 
mentary  Appearance.  Although  it  is  tempting  to 
point  to  the  filamentary  appearance  of  the  emitting 
regions  as  evidence  for  filamentation,  there  are  at 
least  three  alternative  mechanisms  for  explaining 
this  appearance.  First,  since  the  images  are  time- 
integrated,  target  motion  due  to  ablative  accelera¬ 
tion25  during  the  period  of  emission  could  result  in 
an  axial  smearing  in  the  images.  This  is  easily  elim¬ 
inated  by  using  thicker  targets  which  do  not  move. 

Even  with  a  stationary  target,  however,  the 
nj 4  and  nc  density  surfaces  may  move  as  the 
blowoff  plasma  profile  evolves.  To  determine 
whether  this  effect  is  important,  streak  images  of 
the  emitted  light  are  obtained,  using  the  setup 
described  in  Sec.  II.  A  direct  test  for  motion  of  the 
n ,.  surface,  through  streak  photography  of  the  2  w0 
emission,  is  not  possible,  since  the  poor  contrast  of 


the  2  cuq  emission  over  emitted  continuum  prevents 
distinguishing  between  the  two  emissions  without 
two  spatial  dimensions  of  resolution  (one  spatial 
dimension  must  be  sacrificed  for  time  resolution  in 
streak  imaging).  We  do  directly  test,  through  3/2 
<o0  streak  photography,  for  motion  of  the  nj 4  sur¬ 
face,  as  shown  in  Fig.  5.  At  any  given  instant  of 
time,  the  extent  of  the  emitting  region  along  the 
laser  axis  can  be  quite  large,  up  to  250  pm.  The 
position  of  strongest  3/2  w0  emission  is  nearly  sta¬ 
tionary  for  the  entire  emitting  period,  approximately 
2  nsec  in  this  case.  Hence,  the  effect  of  nj 4  surface 
mofiin  on  filament  lengths  is  negligible.  Since  the 
targgt  motion  is  also  negligible,  we  feel  that  the 
motion  of  the  nc  surface,  located  between  the  two, 
is  likely  to  be  unimportant.  A  temporal  modulation 
(approximately  200  ps  period),  as  observed  in  the 
emission  of  Fig.  5,  is  not  uncommon,  and  is  not 
believed  to  be  present  in  the  incident  laser  pulse; 
we  haVe  no  definite  explanation  for  this  modula¬ 
tion,  which  is  reminiscent  of  earlier  results  at 
KMS26'-and  Rutherford27  using  shorter  laser  pulses 
at  higher  incident  intensity. 

The  third  mechanism  which  can  spuriously 
cause  the  filamentary  appearance  of  the  visible 
images  is  refraction.  Refraction  of  the  emitted  light 
as  it  exits  the  blowoff  plasma  can  lead  to  lengthen¬ 
ing  of  the  observed  image,  since  rays  emitted  from 
a  point  into  different  angles  are  bent  differently  and 
appear  to  originate  from  different  sources.  This  is 
illustrated  in  Fig.  6A.  Refracted  ray  paths  are  cal¬ 
culated  for  point  sources  at  nj 4  (for  3/2  o>0)  and  «, 
(for  2  (oq)  in  a  plane-stratified  plasma  with  an 
exponential  density  profile.  For  L„  =  100-125  /zm 
at  nj 4,  the  observed  length  of  3/2  o>0  structures 
could  be  explained  solely  by  refraction.  The 
double-lobe  pattern  of  Fig.  3B  may  also  be  a  signa¬ 
ture  of  refraction.  For  a  certain  range  of  L„,  a  sin¬ 
gle  source  may  emit  two  isolated  sets  of  rays  (one 
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Fig.  S  —  Streak  photograph  of  3/2  w0  emission  from  70  /xm 
thick  CH  target  at  /0  7.2  x  10IJ  W/cmJ.  Filaments  are  ob¬ 

served  in  2-D  image  of  3/2  w0  on  same  shot.  Image  of  target  is 
orthogonal  to  streak  camera  slit,  so  the  one  dimension  of  spatial 
resolution,  labeled  r  in  the  figure,  is  approximately  the  axis  of 
the  incident  beam.  Laser  is  incident  from  right,  and  target  posi¬ 
tion  is  approximately  at  r  ■=  0.  Time  increases  upwards  and 
/  =  0  in  the  figure  is  arbitrary,  though  this  emission  is  occurring 
during  the  peak  intensity  of  the  incident  pulse. 
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POSSIBLE  REFRACTIVE  EFFECTS 
ON  IMAGES  OF  VISIBLE  EMISSION 


(A)  POINT  SOURCE  MAY  APPEAR  FILAMENTARY 

SOURCE:  WMtiM  IMA6E:  W/Z/w/m 


♦ 


IB)  SINGLE  EMITTING  SURFACE  MAY  APPEAR 
OOUBLE-LOBED 


SOURCE: 


IMAGE. 


Fig.  6  —  An  illustration  of  the  possible  difference,  due  to  refraction,  between  fila¬ 
ments  and  filamentary-appearing  structures.  The  hatched  rectangle  is  the  target,  the 
horizontal  lines  below  the  target  are  isodensity  contours,  and  the  laser  is  incident 
from  below  for  all  of  the  illustrations.  (A)  On  the  left,  we  show  a  point  source  of 
2a»0  on  the  nc  contour  or  of  3/2  ai0  on  the  nj 4  contour.  On  the  right,  the  paths  of 
rays  of  harmonic  emission  which  has  been  emitted  into  a  range  of  angles  are  shown. 
The  bending  within  the  plasma  of  rays  that  are  collected  by  the  // 2  lens  leads  to  a 
locus  of  apparent  source  points  on  the  image  which  resemble  a  filament.  (B)  Here, 
the  effect  of  refraction  on  emission  from  a  complete  surface,  as  shown  on  the  left,  is 
examined.  For  a  certain  range  of  density  gradient  scalelengths,  one  can  obtain  the 
family  of  ray  paths  numbered  1-4.  Ray  I,  emitted  near  90°  to  the  incident  laser 
direction,  is  bent  somewhat  but  still  collected  by  the  // 2  lens.  Rays  such  as  2,  emit¬ 
ted  slightly  closer  to  the  target  normal,  are  not  refracted  enough  to  be  collected.  At 
still  steeper  angle  of  emission,  ray  3  penetrates  to  higher  density,  is  refracted  more 
strongly,  and  is  collected  by  the  lens.  Finally,  ray  4  and  other  rays  emitted  near  the 
target  normal  are  refracted  too  strongly  to  be  collected.  The  apparent  source,  due  to 
rays  near  I  and  3,  is  two  lobes  of  emission  separated  by  an  apparently  nonemitting 
region.  This  may  explain  the  appearance  of  Fig.  3B. 
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near  90°  and  one  propagating  toward  critical  den¬ 
sity,  as  shown  in  Fig.  6B)  which  bend  into  the  f/2 
collection  cone.  Refraction  does  not,  however, 
explain  the  difference  at  higher  /0  between  images 
with  the  two  polarizations  (Fig.  3B  vs  3C).  One 
possibility  is  that  two  production  mechanisms 
operate  in  different  regions  of  the  plasma,  with  each 
responsible  for  one  of  the  observed  polarizations. 
Lengths  of  2  <o0  images  could  also  be  explained  by 
refraction,  though  this  requires  an  L„  >  200  /im, 
much  larger  than  expected  at  nc.  Additionally, 
when  the  2  u>0  surface  appears,  it  has  a  much 
smaller  axial  extent  than  the  filaments,  implying 
that  refraction  is  not  the  sole  cause  of  2  w0  filament 
length. 

Inference  of  Locally  Lower  Density.  It  appears, 
then,  that  2  a»0  filament  lengths  are  likely  to  be 
real,  though  refraction  may  be  enhancing  the 
lengths  observed  for  3/2  w0  structures.  If  the  2  w0 
filaments  are  real,  lower  density  may  be  inferred  in 
the  channels;  this  also  applies  for  3/2  w0  if  filament 
lengths  are  not  due  to  refraction.  The  appearance 
of  the  3/2  aiQ  and  2  to0  surfaces  on  the  downstream 
side  of  the  filaments  implies  that  the  nj 4  and  nc 
surfaces  are  perturbed  locally  toward  higher  density 
regions,  as  by  a  low  density  channel  (see  Fig.  7B). 
A  high  density  channel,  as  in  Fig.  7A,  would  per¬ 
turb  the  surfaces  toward  lower  density,  and  they 
would  be  seen  upstream  of  the  filaments.  There¬ 
fore,  visible  observations  do  imply  higher  intensity 
and  lower  density,  as  for  filamentation,  in  regions 
more  localized  than  laser  hot  spots.  That  this  only 
occurs  at  higher  /0  would  also  be  consistent  with 
filamentation,  but  could  be  due  to  the  threshold  for 
production  of  the  harmonic  emission  itself. 


(Al  HIGHER  DENSITY  CHANNELS 


C.  Calculated  Growth  of  Filamentation 

Optical  filamentation  in  plasmas  can  generally 
arise  from  both  thermal3'6  and  ponderomotive3-7'10 
mechanisms.  Under  the  experimental  conditions 
described  here,  these  components  can  have  compar¬ 
able  magnitudes,  although  the  thermal  part  will 
tend  to  predominate  in  higher  Z  targets.  The  fol¬ 
lowing  calculations  will  therefore  include  both 
effects,  and  will  examine  the  applicability  of  the 
theory310  to  the  present  experiments. 

The  calculations  will  proceed  from  three  main 
assumptions:  (a)  the  filamentation  is  quasi-steady- 
state,  (b)  the  plasma  parameters  vary  slowly  in  dis¬ 
tances  comparable  to  filamentation  growth  lengths, 
and  (c)  the  electron  density  variations  n  (r)  -  n0 
induced  by  the  beam  inhomogeneity  remain  small. 
Assumption  (a)  is  justified  if  both  the  ion-acoustic 
transit  time  and  thermal  diffusion  time  across  the 
filament  remain  short  in  comparison  to  the  optical 
pulsewidth.  Although  assumption  (b)  is  not  fully 
satisfied  under  the  present  experimental  conditions 
(plasma  scalelengths  beyond  nj 4  are  on  the  order 
of  one  hundred  microns),  the  calculated  growth 
lengths  are  expected  to  be  a  reasonable  first  approx¬ 
imation.  Assumption  (c)  precludes  any  considera¬ 
tion  of  the  saturation  effects  discussed  by  Sodha,  et 
al.,5  and  effectively  limits  the  theory  to  a  perturba¬ 
tion  treatment.  Within  these  constraints,  the  opti¬ 
cal  behavior  of  the  plasma  is  described  by  the 
refractive  index 

n  =  (1  -  n/n(.)v 2  =  e0l/2  +  1/2  n2  Ifl2,  (1) 

where  e0  =  1  —  no/nc  is  the  linear  plasma  dielectric 
constant,  E  is  the  optical  field  amplitude,  and  n2  is 
the  nonlinear  refractive  index 


(B)  LOWER  DENSITY  CHANNELS 

WiW'/W 


Fig.  7  —  Illustration  of  the  method  used  to  infer  that  harmonic  im¬ 
ages  imply  lower  density.  Target  is  hatched  rectangle,  laser  is  in¬ 
cident  from  below,  and  lines  correspond  to  isodensity  contours.  Solid 
line  corresponds  to  emitting  surface  <n,  for  2m0  or  nj 4  for  3/2  <u0). 
(A)  The  effect  of  high-density  channels  is  to  locally  perturb  isodensi¬ 
ty  surfaces  toward  lower  density  regions.  (B)  Conversely,  lower  den¬ 
sity  channels  perturb  isodensity  surfaces  inward,  toward  the  target. 
This  corresponds  to  the  relative  orientation  of  filaments  and  surfaces 
experimentally  observed  in  the  images. 
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«2  =  2 dh/d  |£|2  =  —  (!///,. «0i/2)  dn/d  \E\ 2.  (2) 


Filamentation  in  a  nonlinear  medium  de¬ 
scribed  by  Eq.  (1)  has  been  treated  theoretically  by 
many  authors.28,6,7  The  basic  approach  is  to  calcu¬ 
late  the  longitudinal  growth  rate  a  (kL)  of  a  small 
transverse  (x)  periodic  perturbation 

£,  Or.  z)  =  E+  (z)  Av  +  £  (z)  <T'V. 

superimposed  on  a  large  plane  wave  E0(z).  Substi¬ 
tuting  the  total  amplitude  £  =  £0  +  Ex  into  the 
optical  wave  equation,  one  obtains  the  general 
expression28 

a  (fc|)=  | (A:x2  w2|£0I2/2€({/2)  -  (/:i2/2/c0)2j  ,  (3) 
where  k0  =  «d/2  w<>/c. 

The  expression  for  n2  must  now  be  calculated 
from  the  plasma  momentum  and  energy  conserva¬ 
tion  equations,  including  ponderomotive  force  and 
collisional  absorption.  Expressing  the  electron  den¬ 
sity  w,  blowoff  velocity  «,  and  electron  (ion)  tem¬ 
perature  Te  (Tj)  in  terms  of  their  quasi-equilibrium 
values,  plus  the  periodic  transverse  perturbations, 

n  =  n0  +  «i  (x,z),  u  =  Z70  +  U\  (x,z) 

Te  =  Tm  4-  T |  (x,z),  Tf  =»  T0  <  Teo, 

one  obtains  the  linearized  steady  state  momentum 
equation 

2 (hq/'Z)  m,  w0  •  V  w,  = 

-  (T„  +  Tq/Z)  V  «,  -  WoV  Ty 

-  (no/16  7T  nc)V  |£|,2.  (4) 

where  |£|2  =  £o  £(  4-  £o£[-  The  steady  state 
equation  is  justified  for  the  conditions  of  interest 
( 7*0  —  eV  and  kL  =  2i T/kL  =  50 /am),  where 
the  acoustic  velocity  cs  is  of  order  2  x  107  cm/sec, 
and  the  transit  time  is  of  order  (kL  c*)-1  =  40  psec. 
In  the  region  near  the  center  of  the  beam  (where 
most  of  the  filamentation  is  observed),  the  blowoff 
is  nearly  longitudinal;  hence,  one  can  take  u0  ■ 
V  M,  =  0.  Equations  (4)  and  (2)  then  yield  the 
result 

w  2  =  n2P  +  n2T,  (5) 

where  n2P  and  n2T  are,  respectively,  the  pondero¬ 
motive  and  thermal  contributions 


"2P  = 

(6a) 

16tt  (1  +  \/Z)  e01/2 

n,  T0  ’ 

n2T  = 

dT\ 

(6b) 

(1  +  1  /Z)  e0l/2  T0 

d\E\r 

and  we  have  approximated  Teo  —  T0. 


The  expression  for  dTx/d\E\2  can  now  be 
found  from  the  linearized  electron  heat  equation 

V  q\  +  «0  T0  V  •  w, 

+  (2 melm,)ve  •  (3/2)  n0  7j 

=  «oI£|2/8  7twc,  (7) 

where  lj\  is  the  thermal  flux  perturbation,  and 

vt  =  5.90  x  10- 24  In  A  ZnJTP  (cgs)  (8) 

is  the  electron-ion  collision  rate.29  The  convection 
term,  which  is  proportional  to  u0  •  V  Tx,  has  been 
neglected  for  a  reason  similar  to  that  given  in  the 
preceding  paragraph  for  neglecting  m0  ’  V 

Following  earlier  treatments,5,6  we  now 
assume  that  qx  can  be  described  classically;30  i.e., 

qx  =  -KV  Tx,  (9a) 

K  =  5.98  x 

1050  r05/2  <t>  (Z)/Z  InA ,  (cgs)  (9b) 

where  =  (Z  +  .24)/  ( 1  +  .24Z).  However, 

according  to  Bell,  Evans  and  Nicholas,31  such  a 
description  will  be  adequate  only  if  |V  Tx\  < 
T(J30\e,  where  \e  =  (TJ me)V2  r"1  is  the  electron 
mean  free  path.  For  T0  —  400  eV,  n0  =  nj 4,  and 
=  50  jam,  one  finds  \e  =  1 1 .5/2  micron;  hence, 
the  classical  description  requires  I Tx\/ T0  <  Z/ 40, 
which  is  consistent  only  with  a  perturbation  treat¬ 
ment. 


One  can  now  show  that  the  thermal  diffusion 
time,  along  with  the  compression  and  electron-ion 
relaxation  terms  of  Eq.  (7),  are  negligible  under  the 
above  considerations.  The  diffusion  time  A  /  can  be 
estimated  from  the  relation  (ve  A  /)l/2X(,  =  kfl, 
giving  the  result  A  t  =  (2.5)Z  psec.  The 

compression/conduction  ratio  can  be  estimated 
from  (9a,b): 


[/7p  7~o  V  •  K)  I 
\K  V2  Tx\ 


X 1/2 
me 


(kL  A,,)-1  «  1. 


Finally,  the  relaxation/conduction  ratio  is 


)3  (mj mj)  ve  n0  Tx\ 

MrV2r,| 


—  (ki\(,)-2«  1 
m. 


Under  the  conditions  of  interest  Eq.  (7)  thus 
reduces  to  the  form5  6 

-K  V2  Tx  =  K  kL2  7j  =  ve  «0  |£|,2/8tt  (10) 

which  gives  the  thermal  nonlinear  index 

-  _ _ <no/nc)2i>e 

"2r  87r  (1  +  \/Z)etf2  k2  KT0 
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when  substituted  into  expression  (6b).  The  relative 
magnitude  of  the  two  indices  is 

"2t/”2P  =  2  Wo  "e/ K  .01 6A|2/A  2<£  ( Z ).  (12) 

This  ratio  is  —  1.8  for  CH  (Z  =  3.5)  and  =16  for 
Al  (Z  =  13). 


Expressions  (3),  (5),  (6a)  and  (11)  can  now 
be  combined  to  give  the  total  longitudinal  growth 
rate 

a  (Aj.)=  |(aPMA;i2/Ao)'f«TM_ (Ai2/2A0)2j  (13) 

where 


k0n2P  Ifpl2  W/;t.)  k0l 

2ed/2  “  4(l  +  l/Z)e03/2c/;(.  T0 

2  92xir, 

(1  +  1/Z)e0  T0  (eV) 

(cm'1) 


(14a) 


(14b) 


would  be  the  maximum  value  of  a  [occurring  when 
kL  =  (2A0aPM)l/2l  if  the  ponderomotive  contribu¬ 
tion  were  acting  alone,7  and 


“TM: 


(cm  ') 


n2T  l£0l2 

1/2 

i,  _ 

(nj nc)2ve  l 

1/2 

2«o/2 

2(l  +  l/Z)ei/2cK  T0 

(15a) 

13.2  In  A 

Z2(n0/nc)iI  (W/cm2) 

1/2 

A  (gm) 

(1  +  l/Z)</>  (Z)«(f2  [Tq  (ev)]5 

(15b) 


would  be  the  maximum  value  (occurring  when 
k2  =  0)  if  the  thermal  component  were  acting 
alone.6  One  can  see  physically  that  as  long  as  steady 
state  conditions  prevail,  the  thermal  contribution 
maximizes  as  kL  —  0  (whole  beam  self-focusing) 
because  thermal  diffusion  tends  to  smear  out  the 
shorter  kl  perturbations.  On  the  other  hand,  tran¬ 
sient  effects  will  tend  to  favor  shorter  AL.  Thus, 
any  incident  laser  hot-spots  are  likely  to  grow  and 
produce  filaments  before  whole  beam  self-focusing 
can  occur. 


In  view  of  the  present  trend  toward  shorter 
wavelength  fusion  drivers,  it  is  worth  noting  that 
the  thermal  contribution  aTM  increases  strongly  as 
A0  decreases.  This  behavior  arises  not  only  from 
the  explicit  A0  dependence,  but  also  from  the  T0'5/2 
term  in  expression  (15b). 

The  total  growth  length  a  1  is  plotted  vs. 
At  =  2 ir/kL  in  Fig.  8  for  the  case  of  CH  and  Al  tar¬ 
gets.  For  Al  targets,  the  growth  is  dominated  by 
the  thermal  self-focusing  mechanism;  for  CH,  the 


Fig.  8  —  Calculated  growth  lengths  a-1  vs  transverse 
wavelength  of  the  perturbation,  assuming  n0“/rr/4, 
T0=  400  eV  and  X  =  1.053  pm.  Solid  lines  pertain  to  CH 
(Z  =  3.5),  and  dashed  lines  pertain  to  Al  (Z  -  13). 


thermal  effects  are  approximately  twice  as  important 
as  ponderomotive  effects.  It  is  evident  that  the 
growth  lengths  at  intensities  of  the  hot  spots  in  the 
incident  beam  are  short  enough  to  account  for  the 
observed  filamentation  and  that,  excepting  for  the 
cutoff  regions  where  a  — *  0,  the  growth  lengths 
vary  slowly  with  Ax.  Since  the  filaments  grow  out 
of  incident  hot  spots  whose  kL  components  lie  in 
this  range,  the  slow  a'1  (A^ )  variation  suggests  that 
their  overall  size  should  not  decrease  as  rapidly  as 
they  would  if  a-1  was  decreasing  as  Ai  decreased. 


D.  Conclusions 

A  comparison  between  the  expected  signa¬ 
tures  of  filamentation  and  our  experimental  obser¬ 
vations  yields  a  substantial  amount  of  evidence  for 
the  presence  of  filamentation.  Localization  is 
implied  by  3/2  wQ  and  2  a>0  images  in  regions  with 
transverse  dimensions  less  than  15  /xm,  as  opposed 
to  the  50-80  fi  m  dimensions  of  incident  hot  spots. 
We  infer  lower  density  locally  from  2  w0  and 
perhaps  3/2  cu0  observations,  and  some  evidence  of 
increased  intensity  modulation  over  that  present  in 
the  incident  beam  is  also  found.  All  of  these 
observations  apply  only  at  higher  incident  intensity. 
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In  addition,  calculations  indicate  that  thermal  self- 
focusing  of  hot  spots  in  the  incident  beam  is 
expected.  This  large  amount  of  circumstantial  evi¬ 
dence  leads  us  to  conclude  that  filamentation  is 
probably  occurring.  One  bright  note  for  laser 
fusion  is  that  the  spatial  scale  of  the  intensity  varia¬ 
tions  caused  in  the  absorption  region  is  very  fine, 
much  smaller  than  the  expected  distance  between 
absorption  and  ablation  regions.  Consequently, 
lateral  energy  transport  between  these  two  regions 
may  well  reduce  ablation  pressure  variations  to 
acceptable  levels.1  There  may  still  be  deleterious 
effects  of  plasma  instabilities  driven  by  the  higher 
intensity  in  the  filaments. 

IV.  OBSERVATION  OF  SPICULES 

As  described  in  Sec.  II,  side-on  images  of  x- 
ray  emission  from  the  laser-irradiated  planar  targets 
are  obtained.  Approximately  half  of  the  pinhole 
images  for  CH  targets  show  fine-scale  structures 
projecting  well  out  in  front  of  the  target  surface,  as 
shown  in  Fig.  9A.  The  finest  observed  transverse 
dimension  of  such  spicular  features  is  pinhole  reso¬ 
lution  limited  at  12  /cm.  When  observed  from  the 
target  rear  at  an  angle  of  45°  to  the  target  normal, 
as  shown  in  Fig.  9B,  these  spicules  appear  to 
emerge  from  small  intense  spots  near  the  target 
surface,  as  those  mentioned  in  Sec.  III. A.  On  most 
images,  a  much  larger  number  of  spots  than  spi¬ 
cules  is  observed.  From  many  of  these  spots,  short 
straight  filamentary  structures  emanate;  spicules. 


which  project  much  further  from  the  target  surface, 
emanate  from  only  a  few  of  the  spots. 

Another  property  of  these  structures  is 
discovered  by  using  thin  targets  which  are  ablatively 
accelerated25  to  speeds  that  allow  significant  target 
displacement  during  the  x-ray  emission  (which  has 
the  laser-pulse  duration25).  Some  x-rays  emitted 
from  the  ablating  front  surface  of  the  target  are 
then  observed  as  the  irradiated  section  translates 
through  positions  behind  the  original  target  location 
(see  Fig.  90.  In  this  image,  one  observes  fine 
scale  striations  along  the  direction  of  target  motion; 
some  of  these  connect  to  spicules  o'  -ved  in  front 
of  the  original  target  position. 

The  angle  at  which  a  spicult  rges  from 
the  blowoff  plasma  depends  upon  i  >sition  rela¬ 

tive  to  the  center  of  the  target.  Sp  •  “manating 
from  source  points  in  the  target  near  the 

center  of  the  x-ray  image  propagate  angles  near 
the  target  normal,  while  those  that  emerge  from 
near  the  edge  of  the  target  curve  away  from  the  tar¬ 
get  normal  at  larger  angles. 

The  spicule  orientations  suggest  that  spicules 
might  be  following  the  fluid  flow  of  material  from 
the  target  surface.  To  test  this  hypothesis,  an 
experiment  is  designed  to  determine  the  pattern  of 
material  flow  from  the  target.32  Aluminum  is 
locally  embedded  in  25  /um  diameter  spots  (up  to 
10  Mm  deep)  within  plastic  targets,  as  indicated  in 


(A)  15/wn  CH  FOIL  (B)  7/*m  CH  DISC  (C)  lOjtm  CH  FOIL 
SHOT  80-9073  SHOT  80-9514  SHOT  80-8968 

Fig.  9  —  (A)  Simultaneous  side-on  images  with  pinhole  diameters 
between  8  and  55  jam.  The  laser  is  incident  from  left  with 
/0  =  2.8  x  I012  W/cm2.  (B)  Images  from  45°  to  rear  of  disk  target 
The  laser  is  directed  out  of  the  page  and  45°  to  the  right  with 
/ o  ~  2.8  x  1012  W/cm2.  (C)  A  side  view  of  an  accelerating  target  at 
/0  —  8.3  x  |012  W/cm2.  The  laser  is  incident  from  the  left  and  some 
x-rays  emanate  from  the  accelerating  target  behind  original  foil  posi¬ 
tion.  The  blackening  in  the  lower  right  is  not  part  of  the  image 
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Fig.  10.  Upon  irradiation  of  this  array  of  spots,  the 
flow  of  ablated  aluminum  from  each  spot  can  be 
tracked  in  side-on  x-ray  images  because  of  the 
strong  Al  x-ray  line  emission  above  1  keV.  An 
example  of  data  obtained  in  this  manner  is  shown 
in  Fig.  11.  The  great  resemblance  of  this  pattern  to 
the  spicule  pattern  suggests  that  spicules  are  follow¬ 
ing  the  fluid  flow  of  material  from  the  target. 


Fig.  10  —  Schematic  illustration  of  experimental  apparatus 
used  in  study  of  fluid  flow  of  ablated  material 


Fig  11  —  Example  of  side-on  x-ray  image  obtained  using 
apparatus  of  Fig  10  Laser  is  incident  from  left  with 
ln  =  .12  x  in12  W/cm2  and  10  jam  thick  CH  target  surface  is  to 
right  Flow  of  Al  from  each  spot  is  observed  as  dark  track  (this 
is  a  negative  image)  of  stronger  x-ray  emission  flowing  to  left 
from  the  target  Darkening  in  upper  left  is  part  of  a  separate 
image 


Spicules  may  also  be  observed  in  images  of 
aluminum  targets,  though  different  pinhole  camera 
filtration  must  be  used.  The  intensity  of  imaged 
radiation  with  the  usual  Be  filters  is  much  greater 
with  Al  targets  than  with  CH  targets,  primarily  due 
to  the  presence  of  strong  Al  line  radiation  in  the 
spectral  band  of  the  camera.  The  dominant  lines  in 
the  spectrum  are  found  through  crystal- 
spectrograph  measurements  to  be  just  above  the 
unionized-Al  absorption  edge  at  hr  =  1.56  keV. 
To  preferentially  filter  out  these  lines,  then,  is 
easily  accomplished  by  adding  an  additional  7 /cm 
Al  foil  filter  to  the  pinhole  cameras.  Images 
obtained  in  this  way  do  show  spicules. 

Previously,14  we  speculated  upon  several  pos¬ 
sible  mechanisms  to  explain  the  presence  of  these 
structures,  some  of  which  depended  upon  the  pres¬ 
ence  of  filamentation.  Since  then,  data  has  been 
obtained  which  casts  a  different  light  upon  our  spi¬ 
cule  observations.  By  using  an  x-ray  crystal  spec¬ 
trograph  with  high  spatial  resolution,  it  has  been 
found  that  the  emission  spectra  of  at  least  some,  if 
not  all,  of  the  observed  spicules  are  characteristic  of 
materials  of  higher  atomic  number  Z  than  carbon. 
This  result  invalidates  the  implicit  assumption  in 
the  discussion  of  possible  mechanisms  in  Ref.  14 
that  the  material  in  the  spicule  is  the  same  as  that 
in  the  blowofT  plasma.  There,  spicules  were 
inferred  to  have  higher  density  than  the  surround¬ 
ing  blowoff  plasma  because  they  were  observed  to 
have  higher  emissivity  and  measurements  suggested 
that  they  had  '.he  same  color  temperature.  How¬ 
ever,  the  presence  of  higher  Z  plasma  renders  inap¬ 
propriate  the  application  of  the  color  temperature 
analysis  and  is,  in  itself,  sufficient  explanation  for 
the  observed  higher  emissivity.  The  observation  of 
spicules  can  be  explained,  therefore,  by  the  pres¬ 
ence  of  microscopic  surface  impurities  on  our 
laser-irradiated  foils,  without  need  for  the  possible 
mechanisms  outlined  in  Ref.  14.  These  impurities 
are  collisionally  confined  in  the  spicule,  much  in 
the  same  wav  as  the  intentionally  introduced  tracer 
materials  in  Figs.  10  and  1 1  are  confined. 

V.  CONCLUSION 

Through  imaging  of  visible  harmonic  and 
subharmonic  emissions  from  laser-irradiated  planar 
targets,  evidence  is  found  for  filamentation  of  the 
incident  beam.  The  fine-scale  of  the  intensity  vari¬ 
ation  in  the  absorption  region,  as  pointed  out  in 
Sec.  Ill,  is  much  less  than  the  expected  separation 
between  critical  density  and  the  ablation  surface; 
lateral  energy  transport  may  then  provide  the 
required  ablation  pressure  symmetry  at  the  target 
surface.  Were  whole-beam  self-focusing  to  occur, 
this  conclusion  would  be  modified;  however,  real 
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laser  beams  and  plasmas  have  nonuniformities,  and 
filamentation  seeded  by  these  is  likely  to  break  up 
the  beam  before  whole-beam  effects  can  grow.  The 
most  severe  effects  of  filamentation,  then,  may  be 
that  plasma  instabilities  driven  by  the  high  intensity 
in  the  filaments  can  grow  and  affect  laser-plasma 
coupling  (through  reflective  instabilities)  or  fuel 
preheat  (through  fast  electron  generation).  Spicule 
observations,  which  were  previously  offered  as  pos¬ 
sible  evidence  for  such  an  effect  of  filamentation, 
are  now  attributed  to  the  presence  of  microscopic, 
higher-Z  impurities  on  our  low-Z  target  foils. 
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